Double-dimer pairings and skew Young diagrams 
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Abstract 



£"> , We study the number of tilings of skew Young diagrams by ribbon tiles shaped 

r j | like Dyck paths, in which the tiles are "vertically decreasing". We use these quan- 

tities to compute pairing probabilities in the double-dimer model: Given a planar 
bipartite graph G with special vertices, called nodes, on the outer face, the double- 
dimer model is formed by the superposition of a uniformly random dimer configu- 
ration (perfect matching) of G together with a random dimer configuration of the 
graph formed from G by deleting the nodes. The double-dimer configuration con- 
sists of loops, doubled edges, and chains that start and end at the boundary nodes. 

\^ ' We are interested in how the chains connect the nodes. An interesting special case 

is when the graph is e(Z x N) and the nodes are at evenly spaced locations on the 

z^» ! boundary R as the grid spacing e — >• 0. 

1 Introduction 

Among the combinatorial objects counted by Catalan numbers C n = (2n)!/(n!(n + 1)!) 
are balanced parentheses expressions (BPEs), Dyck paths, and noncrossing pairings (see 
[Sta99l exercise 6.19]). A word of length 2n in the symbols "(" and ")" is said to 
be a balanced parentheses expression if it can be reduced to the empty word by 
successively removing subwords "()". Equivalently, there are n ('s and n )'s and the 
number of Cs minus the number of )'s in any initial segment is nonnegative (see |Sta99[ 
ex. 6.19(r)]). A Dyck path of length 2n is a map h : {0, 1, 2, . . . , 2n} -> {0, 1, . . . } with 
h(0) = h(2n) = and \h(i + 1) — h(i)\ = 1. There is a bijection between Dyck paths and 
BPEs denned by letting h(i) be the number of Cs minus the number of )'s in the prefix 



of length i (see |Sta99l ex. 6.19(i)]). Dyck paths h are also in bijective correspondence 
with noncrossing pairings n, that is, pairings of points {1, 2, . . . , 2n} arranged counter- 
clockwise on a circle, in which no two matched pairs cross. The bijection is defined as 
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follows: the location of each up step is paired with the first location at which the path 
returns to its previous height before the up step (see |Sta99l ex. 6.19(n)]). 

These sets are also in bijection with "confining" subsets S C {1, . . . , 2n}: a subset S 
is confining if it has the same number of odd and even elements, and for any i 6 S c = 
{1, . . . , 2n} \ S, the set S contains strictly more odds less than % than evens less than i. 
(Confining sets are relevant to the double-dimer model, as discussed in the next section, 
and there the reason for this name will become clear.) The bijection is as follows: odd 
elements in S and even elements in S c are replaced by (; even elements in S and odd 
elements in S c are replaced by ) . The argument is left to the reader. 

For any Dyck path h, let Sh Q {1, . . . , 2n} be its associated confining subset and let 
7th be its associated planar pairing. 

For example, when n — 3, these bijections between the confining sets S, BPEs, Dyck 
paths h, and planar pairings it are summarized in the following table. The pairs in the 
pairing n correspond to the horizontal chords underneath the Dyck path which connect 
each up step with its corresponding down step (shown in the diagrams). The table is 
arranged so that the Dyck paths h are in lexicographic order (equivalently the BPE's are 
in lexicographic order). 



confining set S 



BPEP 



Dyck path h 



pairing n 



diagram 



{1,2,3,4,5,6} 
{1,2,3, 6} 
{1, 4,5,6} 
{1, 6} 



{1, 



3,4, 



6} 



)()() 
)(()) 
())() 
00) 

(())) 



0,1,0,1,0,1,0 



0,1,0,1,2,1,0 



0,1,2,1,0,1,0 



0,1,2,1,2,1,0 



0,1,2,3,2,1,0 



113 15 
6 U \2 



12 3 456 



12345 6 
123 45 6 
12 3 456 



12 3 456 



.0 



We define a binary relation <— on BPEs (and therefore on Dyck paths, etc.) as follows. 
We say P\ <r- P 2 if Pi can be obtained from P 2 by taking some of the matched pairs of 
parentheses of P 2 , and reversing each of them. For example 00 <r- (()) by reversing 
the central pair of matched parentheses. However ()(()) <f- ((())). The relation -(— is 
acyclic and reflexive, and its transitive closure is the well-known partial order ^ on Dyck 
paths, where h\ -< h 2 if h\{i) < h 2 (i) for all i. The lexicographic order used in the above 
table is a linear extension of the partial order -<. 

The following lemma helps motivate the binary relation •(— , and we will use it in 
when we study the double-dimer model. 



Lemma 1.1. Let P±, P 2 be BPEs of equal length, Si the confining subset associated to P± 
and tt 2 the pairing associated to P 2 . Then P\ ^— P 2 if and only if tt 2 has no connection 
from Si to S^. 

Proof. Suppose we reverse pairs pi, ■ ■ ■ ,Pk of P 2 to make Pi. Each pair p^ is a pair of tt 2 . 
If Pi = {a, b}, then a and b have opposite parity. At one of these locations h 2 made an up 
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step, and at the other h 2 made a down step, so either a,b e S 2 or else a, b <£ S 2 . Reversing 
this pair, we have a, b <£ Si or a, b G Si, so the pairs Pi, ■ ■ ■ ,Pk of vr 2 do not match S\ to 
Sf. For any pair {a,b} of 7r 2 that is not reversed, either a,b £ S 2 or else a,b ^ S 2 , in 
which case a, & e Si or else a, b ^ Si, so these pairs do not match Si to S£ either. 

Conversely, suppose tc 2 has no connections from Si to Sf. Consider a pair {a, 6} of n 2 
where a < 6 (so a and 6 have opposite parity and a is an up step of h 2 and b is a down 
step of h 2 ). If a is an up step of hi, then since {a, b} does not connect Si to S£, it must 
be that 6 is a down step of h±. Likewise, if a is a down step of hi, then b is an up step of 
hi. Thus if we reverse the parentheses in BPE P 2 for all such pairs {a, b} for which a is 
a down step in hi, then the result is Pi. □ 

1.1 The incidence matrix M and its inverse M _1 

Associated to the binary relation <— is its "incidence" matrix M, with Mp 1 ^p 2 — 8 o . 
It is convenient to order the rows and columns according the lexicographic order on BPEs. 

Since this order is a linear extension of ^, which in turn is the transitive closure of <—, the 
matrix M will be upper triangular when written in this way. For example, when n = 3, 
Mis 



LO© 


lO"tf 


lOCO 


IH)^ 


lO(M 


co*^ 


co<x> 


CO<M 


CO(M 


CO^ 


i-H(N 


rtCSI 


^H^ 


rH!C 


HtO 


- 


^— s 


•"■"N 


<^*N 


- 



{1,2,3,4,5,6} 


000 


1 


1 


1 


1 


1 


{1,2,3, 6} 


0(0) 





1 





1 





{1, 4,5,6} 


(0)0 








1 


1 





{1, 6} 


(00) 











1 


1 


{1, 3,4, 6} 


((())) 














1 



In addition to labeling the rows and columns of M by the BPEs, we have also labeled 
the rows of M according the corresponding confining subsets, and the columns by the 
corresponding noncrossing pairings, since this is how we shall use the matrices in 



2. 



Since the matrix M is upper triangular with l's on the diagonal, it is invertible and 
M~ l is upper triangular with integer entries. The matrix M _1 is analogous to the Mobius 
function of a partial order, except that the binary relation associated with M is not 
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transitive. For n = 3, M 



is 



to 


to 


to 


to 


to 


LO 




LO 






^~ 




^f 




^ 


CO 


co~ 






CO 


CN 


CN 









000 


1 


-1 


-1 


1 


-2 


0(0) 





1 





-1 


1 


(0)0 








1 


-1 


1 


(00) 











1 


-1 


((())) 














1 



1.2 Skew Young diagrams 

We may associate with each Dyck path an integer partition, or equivalently, a Young 
diagram, which is given by the set of boxes which may be placed above the Dyck path. 
Dyck paths of different lengths may be associated to the same Young diagram, but distinct 
Dyck paths of the same length will be associated to distinct Young diagrams. 

The matrices M and M _1 can be interpreted as being indexed by pairs of Dyck paths, 
which in turn correspond to pairs of Young diagrams. If M AiM is non-zero, then [i is larger 
than A as a path, or equivalently fi is smaller than A as a partition (/i C A). Likewise, 
since M~ l is also upper triangular, if M^ 1 is nonzero then /i C A. Each matrix entry can 
be associated with the skew Young diagram X/fi. Different matrix entries can correspond 
to the same skew Young diagram, but as the next two (easy) lemmas show, when this 
happens, the matrix entries are equal. 

Lemma 1.2. Suppose that fii C \ 1 and fi 2 Q A2, and the skew shapes \\/ [i\ and X 2 / 'fJ.2 
are equivalent in the sense that Ai/ [Xi may be translated to obtain X 2 / V 2 . Then M^ liAtl = 

Proof. Suppose M\ ltfll = 1. Then Ai may be obtained from fii by "pushing down" on 
some of /ii's chords. Each such chord of /j,i must lie within a connected component of the 
skew shape Xi/fii- The result is now easy. □ 

Therefore we may define Mx/^ to be M\ yfJi . 

Lemma 1.3. M^ 1 is determined by the translation equivalence class of the skew shape 

x/n. 
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Proof. We prove this by induction on |A///|. We have 

E % M ft >u 

p:pCpCA 



and isolating the p = X term, 



p-./iQpC-X 



which only depends on A///. □ 



Thus the expression Mr, is well-defined. We give in Figure 1 Mw 1 for the first few 
connected skew shapes. The next lemma shows that for disconnected A///, M^ 1 is a 
product of the M _1 's for the connected components of A//i. 

Lemma 1.4. Suppose fj, C A and A//i /ias fc connected components Ai//^i, . . . , \k/^k- 
Then M x/fl = M Al/m x • • ■ x M Afc/Mfc and M"^ = M" 1 ^ x ■ • • x M" 1 ^ . 

Proof. If we can reverse parentheses in /x's BPE to obtain A's BPE, then in the Dyck 
path representation, the chords connecting the Dyck path steps corresponding to the 
parentheses to be reversed will lie within the region A//x. The multiplicative property 
for Mx/p follows. For Mw , the multiplicative property follows from induction and the 
following equation 

E m *m ■ • • M ^/ Pk M P ~Xi ■ ■ ■ M ;,k = *mi,ai • ■ ■ <w* = <w- □ 



Pi:/*lCpiCAi 
Pk-PkQpkQ^k 



M 



\/p 



o <9 <b <^> 

-11 1-1-2 1 -1 -2 2 -3 -4 -2 2-2 3 

Figure 1: The values of M^} for the first few skew shapes A//i. 

The following theorem gives a formula for M^j . The sign is given by the parity of the 
area of the skew shape. The absolute value is the number of certain tilings, as indicated 
in Figures |2] and [3j Each tile is essentially an expanded version of a Dyck path, where 
each point in a Dyck path is replaced with a box, so we call it a Dyck tile. Dyck tiles 
are ribbon tiles in which the start box and end box are at the same height, and no box 
within the tile is below them. A tiling of the skew Young diagram by Dyck tiles is a Dyck 
tiling. We say that one Dyck tile covers another Dyck tile if the first tile has at least 
one box whose center lies straight above the center of a box in the second tile. We say 
that a Dyck tiling is cover-inclusive if for each pair of its tiles, if the first tile covers the 
second tile, then the horizontal extent of the first tile is included as a subset within the 
horizontal extent of the second tile. We shall prove the following theorem: 
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Figure 2: Cover-inclusive Dyck tilings of the first few skew shapes. 



























Figure 3: Cover-inclusive Dyck tilings of a larger skew shape. 
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Theorem 1.5. M x j = (— l)l A /H x \{cover-inclusive Dyck tilings of X/p}\. 

To prove this, we start with a recursive formula for computing M^} . Using the fact 
M~ l M = J, we get 

E m a:Xm = <w 

E M x/ p M p/, = ^,x- 

p:p<ZpCX 

Let us consider the chords of p that start and end within the region of the skew shape 
X/p. Rather than summing over all p's for which p C p C A, because of the M p /^ factor, 
we can instead sum over all p's which may be obtained from p by pushing down on some 
of the chords which start and end within the region of X/p. 

E m a7p = <w- (i) 

p<ZX:p obtained by pushing down chords of p 

This is the "downward recurrence". We get a different recurrence, the "upward recur- 
rence" from MM' 1 = I: 

E M vM7 = <W 
p 

E M ;l = 6 ^- ( 2 ) 

p~Jp,:\ obtained by pushing down chords of p 

Let us rewrite the downward recurrence for X/fi ^ 0: 

X/p, / j X/ (p with S pushed down) ' \ / 

nonempty sets S of chords of p 

Rather than restrict to sets of chords of \x that can be pushed down to obtain a path p 
above A, it is convenient to sum over all possible (nonempty) sets of chords, with the 
understanding that Mw 1 = if p % X. 





(a) (b) (c) (d) (. 

Figure 4: Illustration of pushing chords down and the resulting Dyck tiles: (a) an example 
upper boundary of a skew shape X/p (the upper boundary need not be a Dyck path, but 
can be extended to a Dyck path), (b) the chords of this upper boundary, corresponding 
to pairs of parentheses that can be reversed, (c) the boundary together with the result of 
pushing down the long chord, (d) the upper boundary together with the result of pushing 
down the middle short chord, (e) if multiple chords are pushed down, the outermost chords 
are pushed down first. Pushing down chords is equivalent to laying down Dyck tiles along 
the upper bondary. 
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Every time we push a chord of fi, we can interpret that as laying down a Dyck tile along 



the upper boundary (adjacent to fi) of the skew shape X/fi (see |Figure 4| ). If multiple 
chords are pushed, then multiple Dyck tiles are laid down, where we follow the convention 
that longer chords are laid down first (so the shorter ones are below the longer ones). If 
we expand the recursive formula for M^} , then each term corresponds to a Dyck tiling 
of X/fi. It is convenient to let each Dyck tile have weight — 1. Since each tile contains an 
odd number of boxes, the parity of the —1 factors in a tiling is just the parity of | A//x| . 
We further rewrite the downward recurrence as 

(-l)l^lM- = £ (-l) 1+l %-l) lX/MS)l M^ iS) l (4) 

nonempty sets S of Dyck tiles 
that can be placed along upper edge of fj, 

where (/i 4- S) denotes fi with S pushed down. 

Let us define fx/p to be the formal linear combination of Dyck tilings of the skew shape 
A//i defined recursively by 

* = Y^ (_i)i+|s| (tiles of 3 placed on top of fx/(^S),\ 

'^ Z-e ^ ' \with longer tiles higher up / ' 



Vwith longer tiles higher up 

nonempty sets S of Dyck tiles 
that can be placed along upper edge of fj, 



(5) 



with the base cases f\/^ = 1 if X/fi = and fx/n = if \i % A. Comparing the recurrences 
OH) and (JSj), we see that if each Dyck tiling is replaced with 1, then f x /^ simplifies to 
(—1)1 VcIm~, . In view of this, we see that ITheorem 1.51 is a corollary of lTheorem 1.61 

Theorem 1.6. With /w^ defined in OS]), fx/n is in fact a linear combination of just the 
cover-inclusive Dyck tilings, with a coefficient of 1 for each such tiling: 

cover-inclusive Dyck tilings T of X/fi 



(See Figure b\ for an example.) 



To prove ITheorem l.Q\ we use the following two lemmas. 

Lemma 1.7. Any Dyck tiling of a skew Young diagram X/fj, ^ contains a tile T along 
the upper boundary of X/fi such that (A//i) \T is a skew Young diagram. 

Proof. Let T\ be the tile of the Dyck tiling which contains the left-most square along the 
upper boundary of X/fi. Suppose tile T n borders the upper boundary of X/fi and has no 
tile above the upper-left edge of its leftmost square (e.g., Ti). Either (X/fi) \T n is a skew 
Young diagram, or else we may define T n+1 to be the tile containing the leftmost square 
square that borders the upper boundary of T n . Tile T n+ i borders the upper boundary of 
X/fi, has no tile above the upper-left edge of its leftmost square, and its leftmost square 
is to the right of the leftmost square of T n . Since there are finitely many tiles in the Dyck 
tiling, the lemma follows. □ 
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♦:♦♦ 



+ ^ +<t> 






-^+«^+<^x+ 



+ 



= + 

Figure 5: Example of using the recursive definition of /a/ m . There is a lot of cancelation, 
and the result is the sum of cover-inclusive Dyck tilings. 

Given a Dyck tiling of a skew Young diagram A/// with k tiles, define a valid labeling 
to be an assignment of the numbers 1, . . . , k to the tiles such that for each j < k, tiles 
1, . . . , j form a skew Young diagram with tile j along its upper boundary. (By lLemma 1.71 
such labelings exist.) Given two tiles T\ and T 2 in the Dyck tiling, say that 7\ -< T 2 if 
Ti's label is smaller than T 2 's label in all such labelings. Then -< is a partial order on the 
tiles of the Dyck tiling. 

Lemma 1.8. Suppose a Dyck tiling of a skew Young diagram A/// contains no pair of 
tiles T\ and T 2 for which (1) T 2 is above T\ (2) T\ and T 2 have no tiles between them, 
and (3) the horizontal extent ofT\ is a proper subset of the extent ofTi. Then the tiling 
is cover-inclusive. 

Proof. We prove the contrapositive. If a Dyck tiling is not cover inclusive, then there 

is a pair of tiles A and B, with B above A, for which the horizontal extent of B is not 

a subset of the extent of A. If the set of squares above A and below B 

(shown in gray in the figure) is nonempty, then let C be any tile containing 

a square between A and B. If the horizontal extent of B is not a subset 

of the extent of C, then we may consider instead the pair of tiles C and 

B, and if the extent of B is a subset of the extent of C, then we may consider instead 

the pair of tiles A and C. For this new pair of tiles, the extent of the upper tile is not a 

subset of the extent of the lower tile, and the interval with respect to the partial order -< 

between the new pair of tiles is strictly smaller than the interval w.r.t. -< between A and 

B. Thus by induction we may assume that tiles A and B have no squares between them. 

If the extent of A is a (proper) subset of the extent of B, then we may take 
T\ = A and T 2 = B. Otherwise, exactly one of the endpoints of A lies within the 
extent of B\ suppose without loss of generality it is the left, and let L denote the 
leftmost box of A (shown in gray in the figure). Since B does not cover A, tile A contains 
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the box immediately up and right of L. Since A and B have no tiles between them, the 
box above L must be part of tile B. The right endpoint of tile B lies above A, and since 
B is a Dyck tile, the box to the immediate upper left of L cannot be part of tile B, say 
that it is part of tile D. This box is the right endpoint of tile D, and since D is a Dyck tile 
and B is a Dyck tile, the left endpoint of B is to the left of D's left endpoint (otherwise 
tile D would have a square lower than its endpoints). Thus, the horizontal extent of D is 
a proper subset of the horizontal extent of B, so we take Ti — D and T 2 = B. □ 

Proof of \Theorem l.b\ We prove the theorem by induction on |A//x|. We have equality 
when X/fi = 0. Otherwise, there is some set S of possible Dyck tiles that may be placed 
along the upper boundary of A//i, which correspond to pushing a single chord of /i down. 
If distinct tiles T ly T 2 G S overlap, then one tile is a subset of the other, say T\ C T 2 (this 
is because two chords of \i cannot have interleaved endpoints). 

Let us evaluate /x/ M restricted to tilings with T 2 at the top edge and T\ directly under 
it. Either T\ and T 2 are pushed down in the same step of the recurrence (jSJ) or T\ is 
pushed down after T 2 . 

The subsets S for the first step of the recurrence in which T 2 is present at the top 
may be paired off with one another so that the symmetric difference of each pair is {Ti}. 
(Here we are not assuming that \x \. S lies above A; when it does not lie above A we have 
/a/045) = °0 Let A,A\J {Ti} be such a pair. Comparing /a/(,4(^u{Ti})) with f\/^ A ) 
restricted to tilings with T 2 at the top edge and T\ directly under it, by the induction 
hypothesis they are exactly the same, but they have opposite signs in the recursive formula 
for fx/fi- Thus fx/p has no tilings in which a tile T\ is directly covered by a strictly longer 
tile T 2 . Now using [Lemma 1.8| we conclude that fx/^ contains only cover-inclusive Dyck 
tilings. 

For a given cover-inclusive Dyck tiling T of skew shape fi/X, let us find its coefficient 
in fx/fj,, which we denote [T\fx/n- Let S be the set of tiles along the top edge of T that 
can be pushed down in the first step of the recurrence. (By lLemma \.1\ S ^ 0.) If any 
tile not in S is pushed down, then the result cannot be extended to T, so we have 

[ti /v, = E(- 1 ) 1+|5 ' l m/v(^) = E(- 1 ) 1+|5 ' 1 = !> 

s'cs s'cs 

S'^0 S'^0 

which completes the induction in the proof of lTheorem 1.61 □ 

Let us define fx/nil) to be the polynomial obtained by giving each tile weight q. Then 

Mx}, = A/m(-1)- (6) 

Given the cover-inclusive Dyck tiling characterization, the next several propositions 
are straightforward to verify. Recall the g-analogue notation 



n-l 



n q = 1 + q H \- q' 

%\ q = n q (n - l) q ..A q 
a \ a\„ 



b) q b\ q (a-b)\ q 

THE ELECTRONIC JOURNAL OF COMBINATORICS 18 (2011), #P130 10 



Proposition 1.9. // the lower boundary of the skew shape X/fi is minimal (V-shaped) 
then fx/M) = g |A/ ^- 




|a/mI 



Proof sketch. There is only one cover-inclusive Dyck tiling; in it each tile has size 1. □ 

Proposition 1.10. // A//i is A-shaped, then f\/^(q) is gl A /^l times a q~ 2 -analogue of a 
binomial coefficient as illustrated in the following example: 




|a/mI 



/min(a, d) + min(5, c) 
\ min(a,<i) 



Proof sketch. In any cover-inclusive Dyck tiling, the tiles with size larger than 1 are all 
A-shaped and centered at the peak on the lower boundary, and their sizes decrease (by 
even numbers) when going up. □ 

Proposition 1.11 (First row of M _1 ). If X is the zigzag path, then f\/^(q) is gl A /^l times 
a product of q^ 2 - analogues of heights, as illustrated in the following example: 




=>- q 52 x l g -2 x 2g-2 x 3 g -2 x 3 g -2 x 4 q -2 x 4 q -2 x 5 g -2 x A q -2 x 3 g -2 x 3 g -2 x 2 g -2 x l g -2 

Proof sketch. In any cover-inclusive Dyck tiling of such 

regions, each tile is a zigzag shape which starts at a 

square with even parity (when the lower- left-most square 

has even parity). Each location with a circled number 

h specifies a number between and h — 1, which is the 

number of squares directly below it that get glued to their lower neighbors to be part of 

the same Dyck tile, as shown in the figure. These tower heights are independent of one 

another and determine the tiling. □ 
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Proposition 1.12. If A/// is a width-one strip, then M^} is up to sign given by a nested 
sequence of products and +1 's, with a +1 for each chord and a times for each minimum, 
as indicated in the figure. 




=> ±[((1 + 1) x (1 + 1 + 1) x (1 + 1 + 1) + 1) x (1 + 1 + 1) + 1 + 1] 

Proof sketch. Each +1 term at a chord (or local maximum) corresponds to the endpoints 
of the chord belonging to the same Dyck tile (which extends no lower than the chord). 
Each x at a local minimum corresponds to making independent choices to the left and 
right of the local minimum. □ 

We also observed some formulas (Conjectures[T]and[2]) for which we do not have proofs. 

Conjecture 1. For any row X of M~ x (of order n), the absolute values of the entries add 
up to a divisor of n\. If the lengths of the chords are measured as 1, 2, 3, ... , then the X th 
row-sum is n\J Y\ chord c t x \ c \, as indicated in the following figure. The q-analogue is 



E 



A/m|/2 



f\/M 



l/2\ 



Til, 



n 



, t ~ 1 chord c of X I 1 1 

[A row specifies the lower path and allows any upper path above it.] 




12! 



71280 



12x1x10x1x1x7x1x4x1x2x1x1 
This formula holds whenever n < 8, and presumably in general. An n = 4 example is 
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1 _l o j. * 2 , o 3 , o 4 , 5 1 X (1 + g ) x (1 + q + q 2 ) x (1 + q + q 2 + q 3 
1 + 2q + 3q + 3q + 2q + q — 



lx(l + g)xlxl 

It is possible to prove this conjecture for the first row, corresponding to the zigzag 
path, for which the row sum is n\. This follows from the formula for the entries of the 



first row (Proposition 1.11), together with a very nice bijection between suitably labeled 
Dyck paths and permutations due to de Medicis and Viennot [dM V94] . The generating 
function for the g-analogue of the first row sum is the g-analogue of a formula of Euler: 



J2nlx n 

n=0 



1- 






X 


1- 




X 




1- 


2x 




2x 
1 



3.r 



the g-analogue of which is 

oo 



x- 



n=0 









1 


1- 






i-qX 


1- 




ql q x 




1- 


q2 q x 




q%x 
q 2 3 q x 



Conjecture 2. For any column \i of M l , the absolute values of the entries add up to the 
product of the heights of the chords under \x, as indicated in the figure. The q-analogue is 

53/a/m(?)= II {height of c) q . 

A chord c of X 

[A column specifies the upper path and allows any lower path below it and above the zigzag.] 




^1x2x2x3x3x3x4x4x4x5x5x6 = 1036800 
This formula holds whenever n < 8, and presumably in general. An n = 4 example is 
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1 + 3g + 5g 2 + 5g 3 + 3g 4 + g 5 = 1 x (1 + q) x (1 + q + q 2 ) x (1 + q + q 2 ) 



2 Double-dimer marginals 



We show how the previous results apply to the double-dimer model. 

Let G be a finite bipartite planar graph with edges having positive real weights, and 

N a set of 2n distinguished vertices on its outer face, which for simplicity we assume 

alternate in color. 

The double-dimer model with nodes N is the probabil- 3 Q , Q ^t,,,,^ .2 

ity measure on configurations obtained by superposing a random 

dimer cover (perfect matching) of Q with a random dimer cover 

of Q \ N (see figure at right). Here by random dimer cover we 

mean a dimer cover chosen randomly for the probability measure 

assigning each cover a probability proportional to the product of 

edge weights of its constituent dimers. Let Z and Z N denote the 

corresponding partition sums. • ohmw 6— — ■ d 

4 1 

In such a configuration the nodes are joined in pairs with lines of dimers. In [KW11] 

we showed how to compute the pairing probabilities, that is, the probability that a given 

pairing of the nodes arises in a random configuration. These pairing probabilities are 



I 



1 1 I 



I U I 



Zjj/Z where Z 



i j 



is 



rational functions of the boundary measurements X^j, where X^j 
the weighted sum of dimer covers of Q \ {i,j} and Z is the weighted sum of dimer covers 
of Q. See [KWllj . For example, the probability of the pairing {{1,2}, {3,4}}, which we 
write as \\\, is X^ 2 X ZA /(Xi )2 X^ A + X 1A X 2 $). 

Knowing how to compute the full distribution on the double-dimer pairings, we wish to 
compute its marginals. Suppose for example that we simply wish to know the probability 
that node i is paired with node j when there is a very large number of nodes, or infinitely 
many nodes. This is an example of a marginal probability question, and computing these 
marginal probabilities can be non-trivial even when the full probability distribution is 
known. In this section we see how to compute these marginal probabilities, and at the 
end we study a natural example that has infinitely many nodes. 
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2.1 Contiguous marginals 

We start by computing the probabilities of subpairings where the nodes of the subpairing 
are contiguous. 

We recall the definition of the quantities D$ for a confining subset S of nodes from 
[KW111 §3]: 

Ds = det [{i tJe8 + ujv) x (-i) ^ 1 -^ 2 ^];: 1 ^; 1 > 

and the fact that 

This allows us to determine all contiguous marginals. For example, row \ \ \ | | of M _1 
yields the marginal 

Pi- 



ll 3 I 5 I \ _ -P{1,2,3,4,5,6} ~ -P{1,2,3,6} ~ ^{1,4,5,6} + -°{1,6} ~ 2-D{i, 3 ,4,6} 

,2Ul6 I" ' ' I ~ 




D 

To derive this formula from the Dyck tilings, we would trans- 

/\ /\ /\ 

late 2 1 1 1 1 ^° a Dyck path 123456, which will be the lower path 
of skew shapes for which we find cover-inclusive Dyck tilings. 
Each tiling corresponds to a term in the formula. The sign 
is negative when the area of the skew shape is odd. The in- 
dices are the locations of the odd-up and even-down steps of 
the upper boundary of the shapes. 



2.2 Local non-contiguous marginals /\/\X -^{1,3,4,6} 



Next we consider marginal probabilities where the set of nodes contained within the chords 
of the marginal are not contiguous, for example 2 I 6 I 4 1 16 I ' ' ' • I n this example there are 
three contiguous blocks of nodes contained within chords: 1-6, 11, and 16. We will show 
how to recursively compute these marginals in terms of marginals with fewer contiguous 
blocks. The previous subsection handled the base case of one contiguous block. 

First observe that if a chord connects two non- contiguous blocks, then there is a finite 
number of ways to "fill in" the space between the two connected blocks, and each of these 
filled in marginals have fewer contiguous blocks. In the example above we have 

Prf 1 l 3 l 5 l n I ^ — PrY 1 I 3 I 5 I 11 I 13 I 15 I ^1 _1_ T^T-f- 1 I 3 I 5 I 11 I 13 I 15 I \ 

rI I 2 I 6 I 4 I 16 I ' ' ' J ~ r L \ 2 I 6 I 4 I 16 I 12 I 14 I* ' ' J "+" rI I 2 I 6 I 4 I 16 I 14 I 12 I' ' ' ) 

Next suppose that each chord connects only nodes within its block. Using the method 
of the previous section, for each contiguous block we associate a formal linear combination 
of subsets of nodes from that block. For example, in the first term above 2 1 1 1 4 1 16 1 12 I il I 
• • • , we have {1, 2, 3, 6}-{l, 6} + {l, 3, 4, 6} for the first block and {11, 16} — {11, 13, 14, 16} 
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for the second block. We may extend the union operator U linearly, and take the union 
over all blocks of these formal linear combinations of subsets of nodes. In this case we get 

({1, 2, 3, 6} - {1, 6} + {1, 3, 4, 6}) |J({11, 16} - {11, 13, 14, 16}) = 

{1, 2, 3, 6, 11, 16} - {1, 6, 11, 16} + {1, 3, 4, 6, 11, 16} 
- {1, 2, 3, 6, 11, 13, 14, 16} + {1,6, 11, 13, 14, 16} -{1,3, 4, 6, 11, 13, 14, 16}. 

Let Yls a sS denote this resulting formal linear combination. For any planar pairing tt, 
we have 



Y^ a sM s , n = < 



1 it is consistent with the marginal of interest 

7i connects each node of a block to another node of the same block, 

but not consistent with marginal 
7i connects a node of some block to a node not in a block 

integer otherwise. 



But pairings tt that fall into the fourth case above can be "filled in" as discussed above, 
and the resulting marginals have fewer blocks and can thus recursively be expressed in 
terms of the Ds's. These can then be subtracted from Yls a s^- 

2.3 Nonlocal marginals 

We have seen how to compute local marginal pairing probabilities. If we wish to under- 



stand the scaling limit of double-dimer pairings such as the ones in |Figure 6\ we need to 
be able to compute marginal pairing probabilities where the nodes in question are well 
separated. For example, we'd like to be able to compute the probability that the upper 
left corner is paired with the middle of the left edge while the upper right corner and 
lower right corner are paired with each other. This is a marginal probability on just four 
nodes, but there are many nodes separating these four nodes. Unfortunately we do not 
know how to compute these non-local marginals in general, even for four nodes out of 2n. 
Experiment seems to indicate that the 4- node marginals depend on all the X^j. 

2.4 Evenly spaced nodes 

As an application of these ideas, consider the double-dimer model on eL x eN, and its 
scaling limit on M. x R + , and suppose that there are 2n boundary nodes which alternate in 



color and are at locations x%, . . . ,X2 n (see |Figure 7 ). In the scaling limit Xij = l/\xi — Xj\, 



(see [KcnOO]), and we have the remarkable formula 

Ds _ TT I _ i(_l)l+i+3 

~n~ ~ 11 ' Xj ' Xi ' 

|KWllt Lemma 5.2]. In |KWllj we proved this formula for finite numbers of nodes on the 
real line bounding the upper half plane, but it also holds for a finite number of nodes on a 
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Figure 6: On the left are the chains in a random double-dimer configuration in a box- 
shaped region with many nodes all along its boundary. The loops and doubled edges of 
the double-dimer configuration are not shown here. In the figure on the right, the box 
was conformally mapped to the disk, and the chains connecting points on the circle are 
replaced with arcs. We expect that in the scaling limit, such random chord diagrams will 
be Mobius invariant. 



circle bounding a disk. We can use the right-hand side to define Ds/D^ for non-balanced 
sets S of the nodes. Using 



D 



W 



Do 



n 



tXj rt tXj r> 



(_l)l+*+J 



we can rewrite 



Ei 



D 11 l^i •'./ 



ies ® ijeS 

i<j 



\2(-l) i+ i 



Now suppose that the domain is the unit disk and there are 2n nodes, one at each of 
the 2n th roots of unity. Let ( = e lm l n and u> = e 7 ™/™ . 



IK*-? 

3=1 



z n -l = (z- l)(z n - L + z n ~ 2 + --- + Z + 1 



Thus 



Du } nu\"-c j 



D * nsii-c 



3\ n 
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Sis 




SLlligl 



Figure 7: A random double-dimer configuration with nodes of alternating color at evenly 
spaced locations. We are interested in the limit where the domain converges to the upper 
half plane, the lattice spacing tends to 0, and there is a node at each integer. 



We also have 



Hence 



Jb i Jb ij £i 



sin 



n 



sin 



D a n\ s \ 



ijeS 

1 j — i even 



2n 

2n 



n 

i,jes 

Kj 
j — i odd 



sin 



In 



For example, if j — i is odd then 

D{i, j} 4 1 



1 1 



Do 



n 2 \xi — Xj\ 2 n 2 s i n 2 ?El? H 



-»• 



4/tt 2 



»-J 



in the n — >■ oo limit. More generally, for any finite balanced set 5, in the n — > oo limit 



n (*-i) s 






\S\ »<?' 
j — i even 



n 

i,jeS 

Kj 
j — i odd 



*-J 



For example, the probability that 1 pairs with 2 and 3 pairs with 4 is 

'in / / \ i 

0.246979 



p r r. . . 1 1 1 3 I 



£{1,2,3,4} £{1,4} /2 V 16 

£* D„ U/ 9 



nj 9 
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3 Pairings in groves 

Many results for dimers have analogues for trees and vice versa 
(see e.g., |Tem81l IKPWOOl IKS03J ). The analogue of the double- 
dimer pairings of the previous section are "grove pairings." Given 
a graph (not necessarily bipartite) with a specified subset N of 
the vertices called nodes, a grove is a forest such that each tree 
contains at least one of the nodes (see figure at right). A grove 
pairing is a grove in which each tree contains exactly two such 
nodes. If all such groves have probability proportional to the 
product of their edge weights, we are interested in the probability 
that a random grove has a specified pairing type. 

An important formula for grove pairings, which is due to Curtis, Ingerman, and 
Morrow [CIM98] (see also Fomin |Fom01j ). expresses grove pairing probabilities in terms 
of determinants of a submatrix of the "(electrical) response matrix". If node i is held 
at 1 volt while the other nodes are held at volts, then Lij denotes the current flowing 
out of the network through node j. Though not obvious from this definition, L 
(see e.g., |CdV98j ). If there are 2n nodes {v\, . . . 
Morrow formula is equivalent to 




'-.; 



J .i-i 



Wi, . . . , u> n }, the Curtis- Ingerman- 



det [Lij] 



l=Vl,...,V n 

j=w 1 ,...,w n 



E 

ctG5„ 



•lr- 



z 



"l.^(l)|->«A(n) 



7 ' 

Zj vi\---\v n \wi\---\w„ 



where Z n is the weighted sum of groves where the nodes are connected according to the 
partition n. This formula holds for any graph, whether or not it is planar. In the case 
of planar graphs, we require the nodes to be a subset of the vertices on the outer face, 
which we number in circular order. In this case, if all the v $'s are contiguous, then there 
is only one pairing between the v's and w's for which the term on the right-hand-side is 
nonzero, so for this pairing there is a simple determinant formula. 

When the graph is planar and the nodes are on the outer face, there are C n = ( 2 ^) -^ 
noncrossing pairings n of the nodes, and ( ™)/2 formulas relating the Z„-'s coming from 
the CIM formula, so Dubedat conjectured that all of the Z 7r /Z 1 |...| 2n 's are determined by 
the CIM determinants |Dub06j . Here we show how to adapt our double-dimer results from 
the previous section to prove this conjecture. Essentially the same formulas hold as for 
the double-dimer pairings, except that there are some extra minus signs in the formulas 
for grove pairings. 

Let S* = {v i, . . . , v n } C {1, . . . , 2n}. Then the CIM determinant enumerates (with 
sign) grove pairings compatible with S*, where a pairing is compatible with S* iff every 
chord connects an element of S* to an element not in S*. Given S*, we can define S by 

S = (odds of S*) U (evens not in S*). 

Then the pairing is compatible with S iff no chord connects an element of S to an element 
not in S, which is precisely the same notion of compatibility that we used in the previous 
section on double-dimer pairings. 
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We can rewrite the previous equation as 



det [ L iA%% 



V A|.-|2n 



since whenever Z n ^ 0, we have Ms )7T = 1. Here <j(tt, S*) is the permutation mapping S* 
to the complement of S*, when these sets are listed in sorted order. Conveniently, this 
sign (— l) ^ 71 "' 5 "") does not depend on S*, so we can write it as (— 1)°": 

Lemma 3.1. Let n be a planar pairing. Suppose S* is a set of nodes that is compatible 
with 7r in the sense that each chord connects a node of S* to a node of (S*) c . If the nodes 
of S* are arranged in sorted order and likewise for (S*) c , then the permutation defined by 
the pairing has sign which is determined by is alone and is independent of S*. 

Proof. It suffices to compare the signs for £* and S% which differ at one chord {a, b} of 
the pairing, where say a < b and a & SI and b G S%. When a G S± is replaced with b, it 
is moved to the right a number of times for the list to remain in sorted order, and this 
number equals the number of chords {c, d} where a < c < d < b. Likewise when b G S* is 
replaced with a, it is moved to the left the same number of times for the list to remain in 
sorted order. Thus the parity of the permutation connecting SI to {Sl) c changes an even 
number of times when it is transformed to the permutation connecting S^ to {S^Y- □ 

Therefore we can solve for the Z„-'s in the same way that we solved for the Pr|Vr]'s in 
the previous section: 

Theorem 3.2. If ir is a noncrossing pairing of {1, . . . , 2n}, for groves in a planar graph 
we have 



:-i) 7 



ZiH 



2/i 



zZ M ^det[L id ]%%. 



113 15 



For example, when there are 6 nodes, the partition function for the pairing \ \ A (j 
given by 



is 



Zl\3\5 
2UI6 



-Zl|2|3|4| 



5 6 



Li : 2 L\, 
L-i L 



^3,2 
-^5,2 L 



3,-1 

5,4 



u 


1,6 




-^1,2 


-^1,5 


L 


1,6 










£3,6 


— 


£3,2 


-^3,5 £3,6 






£5,6 




-^4,2 


^4,5 -^4,6 








^1,3 


Ll t 4 


Li,6 




-^1,3 


-^1,5 -^1,6 




-^1,4 -^1,5 -^1,6 




-^2,3 


-^2,4 


-^2,6 


+ 


-^2,3 


-^2,5 -^2,6 


-2 


-^2,4 -^2,5 -^2,6 




L 


5,3 


-^5,4 


-^5,6 




L 


1,3 


£4,5 -^4,6 




-^3,4 L 3i5 L 3i6 
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To derive this formula from the Dyck tilings, as with the 
double-dimer pairings, we translate \ \ \ | | to a Dyck path 
123456, which will be the lower path of skew shapes for which 
we find cover- inclusive Dyck tilings (see figure to right). Each 
tiling corresponds to a term in the formula. The sign is neg- 
ative when the area of the skew shape is odd. The rows are 
indexed by the locations of the up steps of the upper bound- 
ary, and the columns are indexed by the down steps of the 
upper boundary. In addition there is a sign on the left-hand- 
side which is the sign of the permutation associated with the 
pairing. 

If we cyclically rotate the indices, then the pairing \\\\\ 

translates to the Dyck path 23456 1 , which allows us to write 
an equivalent formula with just two determinants, as shown 
in the figure to the right. 

Note that we only know how to compute the grove partition function Z n for a complete 
pairing 7r; unlike the situation for double-dimer pairings, we do not know how to compute 
marginal probabilities. 





+ det [L itj ] 
-det [Li j] 1 . 



«=2,3,5 

J=4,6,l 

2,3,4 

:5,6,1 



References 



[CdV98] 

[CIM98] 

[dMV94] 

[Dub06] 
[FomOl] 

[KenOO] 
[KPWOO] 



Spectres de Graphes, volume 4 of Cours Specialises 
Societe Mathematique de France, Paris, 1998. 



Circular planar graphs 
283(1-3):115-150, 1998. 



Yves Colin de Verdiere. 
[Specialized Courses]. 
|MR1652692 (99k:05108)j 

E. B. Curtis, D. Ingerman, and J. A. Morrow, 
and resistor networks. Linear Algebra Appi, 
|MR1657214 (99k:05096)| 

Anne de Medicis and Xavier G. Viennot. Moments des g-polynomes de Laguerre 
et la bijection dc Foata-Zeilberger. Adv. in Appi. Math.. 15(3):262-304. 1994. 
|MR1288802 (95i:05007)j 

Julien Dubedat. Euler integrals for commuting SLEs. J. Stat. Phys., 
123(6):1183-1218, 2006, |arXiv:math.PR/0507276| |MR2253875 (2007g:82027)1 

Sergey Fomin. Loop-erased walks and total positivity. Trans. Araer. 
Math. Soc., 353(9):3563-3583 (electronic), 2001, |arXiv:math.CO7 0004083. 



MR1837248 (2002f:15030) 



Richard Kenyon. Conformal invariance of domino tiling. Ann. Probab., 
28(2):759-795, 2000, |arXiv:math-ph/99i0002J |MR1782431 (2002e:52022) 

Richard W. Kenyon, James G. Propp, and David B. Wilson. Trees and 
matchings. Electron. J. Combin., 7:Research Paper 25, 34 pp., 2000, 



|arXiv:math. CO/9903025 . MR 1756162 (2001a:05 123) 



THE ELECTRONIC JOURNAL OF COMBINATORICS 18 (2011), #P130 



21 



[KS04] 



[KW11] 



Richard W. 
J. Combin. 



Kenyon and Scott Sheffield. 
Theory Ser. B, 92 (2): 295-3 17, 



MR2099145 (2005k:60031) 



Dimers, tilings and trees. 
2004, |arXiv:math/0310195| 



Boundary partitions in 
363(3):1325-1364, 2011, 



Richard W. Kenyon and David B. Wilson. 
trees and dimers. Trans. Amer. Math. Soc, 
|arXiv:math.PR/0608422[ 

[Sta99] Richard P. Stanley. Enumerative Combinatorics, Volume 2. Cambridge Studies 
in Advanced Mathematics #62. Cambridge University Press, Cambridge, 1999. 
With an appendix by Sergey Fomin. MR1 676282 (2000k:05026)1 

[Tem81] H. N. V. Temperley, editor. Combinatorics, volume 52 of London Mathematical 
Society Lecture Note Series, Cambridge, 1981. Cambridge University Press. 
MR633 646 (82i:05001)| 



THE ELECTRONIC JOURNAL OF COMBINATORICS 18 (2011), #P130 



22 



